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Abstract 


A theoretical and experimental investigation of the creep 
behavior Of initially flat circular plates is described in 
which particular attention is paid to the effect on this 
behavior of geometry changes induced during deformation. A 
method of modeling the constitutive behavior of a material 
undergoing creep deformation in a multiaxial state of stress 
is proposed. An energy technique for solving creep problems 
using this constitutive model is described and solutions are 
obtained for the cases of clamped and simply supported circular 
plates under a uniform, time independent, lateral load. The 
technique retains the effects of both bending moments and in- 
plane forces in the deformation of the plate. A series of 
experiments with lead plates is» described and compared with 
the theoretical predictions obtained in this study and with 
those of previous investigators. Both the theoretical and 
experimental results show that the effect of finite deformation 
is of primary importance in the creep behavior of plates and 
must be taken into consideration even when the maximum 
deflection is well below the plate thickness. 
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INTRODUCTION 


Creep is a term commonly applied to a material deformation 
process in which the strains are a function of time as well as 
stress. The phenonmenon occurs in metals when they are stressed 
at high temperatures, and can result in appreciable plastic 
deformation even for stresses well below the yield stress. As 
machinery operating temperatures have risen steadily over the 
past years in search of higher efficiencies the phenomenon of 
creep has been of increasing interest to engineers. When 
designing components which will operate in an environment where 
creep may occur it is obviously necessary for the engineer to 
be able to predict creep deformations in order to produce an 
efficient and workable design. The mechanisms of creep are 
complex and as yet imperfectly understood, but enough is known 
to permit a reasonable analysis of some cases of practical 
importance. It is the purpose of this investigation to examine 
one such case, that of a thin circular plate exposed to a 
uniform lateral pressure, and to egies a means of predicting 
its creep behavior including, in particular, the influence 
upon this behavior of the geometry changes resulting from the 
deformation. It is believed that the general procedure 
developed in the solution of this particular problem can be 
applied to other structural shapes. 

Trouton and Rankine observed and described the phenomenon 
of creep in lead wire in 1904 (1), In 1910, Andrade (2) made 


the first systematic investigation of creep using wires of 





moproximately pure lead and of lead-tin alloy. In 1929, 

Norton (3) published the results of his investigation of creep 
in steel at elevated temperatures and suggested that creep 

rate, for uniaxial stress, could be represented as a function 
of stress raised to a power n. This relation, which has 

become known as Norton's law, has been the basis for much of 
the phenomenological creep analysis in recent years. Other 
functional relations for uniaxial test data have been suggested, 
such as a hyperbolic sine function proposed by Prandtl and 
endorsed by Nadai (4) and an exvonential relation suggested by 
Soderberg (5). The greater mathematical tractability of 
Norton's law has made its use far more common. A generalization 
of Norton's law to the multidimensional stress state was pro- 
posed by Odqvist in 1933 (6). Odqvist used the concept of 
Octahedral stress to relate the 5 eas rate components to the 
deviatoric stresses. At about the same time, Bailey (7) 
proposed a somewhat more complex set of relations for general- 
izing Norton's law to the multiaxial case in a form appropriate 
SOmthne principal directions of stress and strain. Soderberg (5) 
showed how his exponential creep rule could be generalized to 
multidimensional cases through the stress and strain invariants. 
Another approach, based on the work of Drucker (8) and Hopkins 
and Prager (9) in the rigid plastic behavior of circular plates, 
uses the Tresca yield surface and the condition of normality of 
Strain rate vectors to that surface to extend a uniaxial creep 
law to a multiaxial problem. This method was used by Wahl (10) 


for the analysis of rotating discs and by Venkatraman and Hodge 
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(11) for the bending of circular plates subject to creep 
deformation. At the present time there seems to be no firm 
basis, experimental or otherwise, for endorsing or rejecting 

any of the above theories as a generally applicable constitutive 
relation in creep. In this paper another model for the general- 
ization of Norton's law to multidimensional cases is proposed 

in which principal strain rates are assumed to be a function of 
the nth powers of the principal stresses. 

The analysis of the creep behavior of plates has thus far 
been confined to a bending or infinitesimal deflection analysis, 
in which in-plane forces are ignored, or a membrane analysis, 
in which moments are ignored. Odqvist (12) analyzed the bending 
of thin circular plates with clamped and simply supported edges 
uSing his multiaxial relations based on stress invariants. 
Venkatraman and Hodge (11) solved the same bending problems for 
circular plates using the Tresca criterion technique mentioned 
above. Venkatraman and Patel (13) (14) used this same method 
for the bending of annular plates. The Tresca technique was 
extended to the case of orthotropic bending of circular plates 
by Sandaranarayanan (15). A membrane analysis for circular 
plates undergoing creep deformation has been done by Odqvist 
(12) (15), again using his stress invariant generalization of 
Norton's law. Onat and Yuksel (17) have suggested a very 
Simple and ingenious technique for the membrane analysis of 
a thin circular plate in which it is assumed that the plate 


deforms through creep behavior into a spherical shape. The 
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stresses may then be assumed to be everywhere equal and the 
uniaxial stress-strain rate rule is applied directly. There 
has apparently been, as yet, no mathematical analysis of the 
creep behavior of plates which retains both bending moments 
and in-plate forces and which considers the influence of the 
geometry changes induced by the deformation. None of the plate 
analyses cited above were accompanied by any experimental 
verification nor have any creep tests on plates been noted by 
the author in the literature. Jones (18)(19) has shown the 
very important influence of finite deformations in the static 
and dynamic behavior of rigid-plastic plates and it would 
seem likely that creep behavior would exhibit a similar 
devendence. 

In this investigation a series of initially flat lead 
plates of three different thicknesses is subjected to various 
lateral hydrostatic pressures in a test apparatus which permits 
measurement of plate deflection at the center and at regular 
intervals along the radius. By observing deflections at 
various time intervals one is able to determine the deflection 
rate of the plate and the plate's deformation profile. Tensile 
creep tests are also conducted on specimens made from the same 
lead sheets as the circular plates in order to determine basic 
material creep parameters. Experimental procedures, apparatus, 
and results are described in Chapter I. A method for mathe- 
matically analyzing the creep behavior of circular plates 
including the effect of finite deformations is described in 


Chapter II. The method is based upon computing the rate of 
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micernal energy dissipation in the plate and equating it to the ex- 
ternal work rate. It is essentially an adaptation of a technique 
suggested by Timoshenko (20) for computing the large deflection 
behavior of circular elastic plates. The problem under study 
has several inherent nonlinear aspects: nonlinear material 
constitutive relations, nonlinear geometry with finite 
deflection, and nonlinear change in creep parameters with time. 
Despite these multiple nonlinearities, by using the technique 
herein described, the results for a circular plate may be 
presented in the simple form of an algebraic equation and a 

set of graphs. All necessary computer work is incorporated in 
the graphs from which the deformation rates may be computed 
directly, with the aid of a slide rule, for plates made from 
Meyemacerial whose creep characteristics are adequately 
described by Norton's law. In Chapter III the theory is 
compared with the experimental results and with the bending 

and membrane analyses of other investigators. 

Both the theoretical and experimental results demonstrate 
the strong influence of finite deflections on the creep behavior 
Sec ircular plates. This influence is significant even for a 
deflection at the center less than the plate thickness. The 
proposed method of analysis yields results which are in quite 


reasonable agreement with the experimental data. 
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CHAPTER I 
A. PLATE DEFORMATION EXPERIMENTS 
fm DGhoCRIPTION OF APPARATUS 
The creep behavior of a series of circular plates of 
very nearly pure lead was examined in the apparatus shown in 
Figures l and 2. The plates were cut from commercially 
rolled sheets of three different thicknesses commonly known 


as 4 lb., 6 lb., and 8 lb. plate. These thicknesses were: 


Plate Size Thickness 
4 lb 2060" 
6 lb ad 0 oo a 
sale: esiO e 


Local variations within the plates from the above thicknesses 
were no greater than +.001". All circular plate specimens and 
tensile specimens of any one thickness were cut from the same 
lead sheet. A chemical analysis of each sheet is given in 
Appendix G. Lead was chosen for these experiments because of 
the fact that it creeps at room temperature. This fortunate 
property means that costly furnaces, temperature control 
devices, and elaborate instrumentation which would be necessary 
for elevated temperature testing are not required. Although 
lead is not a commonly used structural material we may expect 
that its behavior at room temperature will correspond to the 
behavior of more common structural materials, such as steel, 
at higher temperatures. Moreover, the theoretical analysis 
which is to be presented is valid for any material whose creep 
behavior is governed by the constitutive equation postualted. 


Thus, 1f agreement between theoretical predictions 
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and experiments is achieved for lead plates, then the analysis 
should remain valid for plates which are made from other 
materials whose behavior is adequately represented by the 
constitutive equation. 

The test apparatus consisted of a short section of steel 
pipe with a blank base welded to one end, and a flange and 
ring on the other end which was used to clamp the test 
Specimen. The flange and ring were machined to have an 
inner diameter of exactly five inches and square edges at the 
inner lip. Centering pegs were placed in the flange so that 
the ring could always be positioned precisely above it by 
fitting the pegs into holes drilled into the ring. This was 
necessary as the inner edge of the ring must be exactly above 
the edge of the pressure vessel to insure proper end conditions. 
The flange and ring had serrations machined on the gripping 
Surfaces. These serrations were relieved in both surfaces so 
as not to damage or distort the edge of the plate sample. The 
Sample plate was clamped between the flange and the ring by 
eight high strength steel bolts tightened to a uniform tension 
with a torque wrench. Six dial indicators were positioned 
above the plate with indicator holders attached to the ring. 
These dial indicators were placed so as to record the defor- 
mations at various places on the plate. Positioning of the 
indicators was facilitated by the use of extension rods of 
various lengths for the indicator contact points. All dial 
indicators were graduated to .0001" and were checked for 


accuracy with gage blocks before commencing the experiments. 
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The apparatus was pressurized by a water column. The water 
column consisted of a length of plastic tubing connected to a 
3" diameter plastic cylinder at one end and, through a vaive, 
to the pressure vessel at the other end. The use of a water 
column insured that the hydrostatic pressure in the vessel 
could be both accurately measured and accurately maintained 
over the long period of time required for a creep test. The 
meeender on the upper end of the tube was for the purpose of 
expanding the volume of the upper part of the column so that 
small changes in the volume of the pressure vessel from plate 
deformation could not make an appreciable change in the column 
height. The pressure vessel was also equipped with a vent 
valve and a pressure gage. The pressure gage was used as a 
check to insure the pressure in the vessel was at the level it 
Should be and not blocked in some way. Test pressure data 

was obtained by measuring the water column height not by reading 


the gage. 


2. EXPERIMENTAL PROCEDURE 

The circular plate specimens, having been cut from a 
Single sheet of commercially rolled lead, were flattened with 
a hydraulic press and then annealed at 200°C for two hours. 
After cooling to room temperature a sample plate would be 
mounted in the test apparatus as shown in Proiren bea. e The 
Chamber of the test apparatus and the pressurizing column 
were filled with water before putting the sample in place. 


The sample plate was marked by pencil with concentric circles 
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from the center at 1/2 inch intervals, and then placed 
carefully on the upper flange of the test apparatus. The 
securing ring was then put into position above the specimen 
and secured with eight bolts tightened to a uniform tension 
with a torque wrench. The dial indicators, having been 
previously attached to the securing ring, were then positioned 
to their desired locations at the center and at intervals 
along the radius. The assumption was made that deformation 
was axisymmetric so that readings would be a function of 
radius only. Observation of the tests revealed no evidence 
that this was invalid. With the water column adjusted to a 
height giving the desired hydrostatic pressure the test is 
ready to begin. Before releasing the pressure, initial 
readings are taken on all dial indicators. These are the flat 
olate, or t=0, readings to which all subsequent readings are 
referred. The test is started by opening the valve between 
the pressure column and the pressure vessel, thus admitting 

a hydrostatic pressure to the back of the plate. Dial 
indicator readings and room temperature are subsequently 
recorded at various time intervals. Plotting these readings 
at any given time aives the deformation profile at that time, 
and plotting the readings at any given radius over all times 
gives the deformation history. The slope of the deformation 
history at any particular time is the deformation rate at that 
time. Of particular interest is the largest deformation which 
occurs at the center. Two tests were made at each pressure 


using a different plate of the same thickness in order to check 


the reproducibility of the results. 
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Pe RESULTS 

The results of the plate deformation experiments are 
summarized in Figures 1.4 through 1.16. Raw data is listed 
in Avpendix F. 

Figures 1.4 through 1.6 are plots of the deflection of 
the center of the plate as a function of time. Two ordinate 
scales are shown, one for the deflection in inches and 
another for deflection as a ratio of the corresponding plate 
thickness. 

Figures 1.7 through 1.9 also show center deflection as a 
function of time, but here all deflections are plotted relative 
to that at t=1l hour in order to show more clearly the relative 
deflection rates for different meee enr es. The time t=1l hour 
was chosen arbitrarily as a reference point. Some other point 
could have served equally well. 

Figures 1.10 through 1.16 show representative deflection 


Ww 
° e s O 
Meeritles for increasing values of —. The plots are of trans- 


h 
verse deflection as a function of radius. Those shown are at 
time t=1l hour. Similar plots for other times such as t=5 hours 
Show no significant differences. It can be seen from these 


plots that the deflection profile is approximately bounded by 


curves of the forms 


hee 2 
vo = Ww (l-=z) dak 
and 
r? 
Ww = we oe) ee 


where we is deflection at the plate center. For tests at the 
Ww 


lower = ratios, the measured profiles correspond quite closely 
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W 
to the former expression, and at higher —° ratios the observed 


h 
profile seems to migrate toward the latter parabolic form. 
This seems rather reasonable as the former expression corresponds 
to the profile given by Timoshenko (21) for the elastic 
deformation of a clamped, hydrostatically loaded, circular 
plate and the latter is the corresponding membrane profile. 
The fact that the plate behaves in this way was of particular 


Significance in suggesting the method of analysis described in 


the next chapter. 


poe LENGLLE CREEP TESTS 
1. DESCRIPTION OF APPARATUS 

Tensile creep tests were necessary for obtaining the 
appropriate creep constants, that is, K and n in Norton's law 
which may be written as: 

é = Ko" 5 
Tensile test specimens were prepared from the same lead sheets 
from which the plate specimens were made. The tensile 
specimens were also flattened with the same hydraulic press 
and annealed for two hours at 200°C in a manner similar to 
the circular plates. The type of specimen used is illustrated 
i FiqQuyew ik? 

The tensile testing apparatus is shown in Figures 1.17 
A, B, C, and EF. It consists essentially of a frame in which 
the weighted tensile specimen is suspended and a cathetometer, 
an optical measuring device, for measuring the distance 


between two marks inscribed on the specimen. There is also 
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a modified hydraulic jack used to facilitate the rather 
important task of smoothly loading the specimen at the 
commencement of a test. The specimens were loaded with 
calibrated iron weights mounted on a pan suspended from the 
lower s~ecimen's grip. 

The frame was constructed of aluminum I beams bolted 
together. The specimen was suspended from an eye mounted 
in a transverse bar centered in the top of the frame. A 
flexible braided wire with an eye soldered to each end was 
used to suspend the specimen. Each of these eyes had its 
inner edge machined to a knife edge to insure that a specimen 
would hang freely. A slotted bar, as shown in Figure 1.17C, 
was used to prevent any rotation of the specimen during the 
test. This was necessary to insure the marks were always 
visible by the cathetometer. The bar was placed loosely on 
its supports and provided only a very gentle, yet sufficient, 
restraint against any tendency for the specimen to rotate. 

The specimen was held between two grips as shown in 
Figure 1.17C. The upper grip was made of aluminum and the 
lower of plastic. The use of plastic was to minimize the 
weight of the lower grip. This assisted in preventing any 
bending of the specimens while mounting them in the frame. 

The initial loading of the specimens was accomplished 
Peeieene aid of a modified hydraulic jack. Prior to the 
commencement of a test the weights and weight pan were 
supported by a small hydraulic jack as shown in Figure 1.17E. 


The jack is adjusted so that the hook on the weight pan is 
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barely out of contact with the ring on the lower grip. When 
the jack's bypass valve is released, two springs draw the 
piston quickly and smoothly down from the weight pan, thus 
loading the specimen. As seen in Figure 1.17B the radial 
slots in the weights ride on a verticle rod. This rod 
prevents the weights from rotatina and twisting the specimen 
but does not inhibit vertical movement of the weights. 
Measurement of tensile strain was accomplished with the 
aid of the optical device shown in Figure 1.17 A, B, and C, 
known aS a cathetometer. It consists of two telescopes 
mounted so that they can be moved vertically as a unit or 
moved apart in the vertical direction. The unit had a vernier 
for measuring a change in the vertical position of the two 
telescopes aS a unit. Each telescope was equipped with a 
movable hairline whose deflection from its centered position 
was measured by a vernier graduated to .0001". By following 
the separation of two scribed marks initially a known distance 
apart on the specimen, one can use this instrument to measure 
creep elongation and have no physical contact with the 
specimen which could possibly affect the normal deformation 


process. 


2. EXPERIMENTAL PROCEDURE 

The tensile specimens, having been prepared from the lead 
Sheet as described in the previous section, were scribed with 
two thin marks one inch apart on the narrow section of the 
sample. The scribe marks were made by clamping a one inch 


gage block over the specimen and drawing a razor blade down 
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either side of the gage block. Great care was necessary in 
attaching the grips to the specimen in order to avoid bending 
the soft lead. A svecially prepared block, machined to the 
proper dimensions, was found useful in this regard for 
supporting the specimen and grips while they were being 
attached. The slotted bar, used to prevent rotation of the 
sample, was placed around the specimen when attaching the 
griovs. This assemblage is then suspended in the frame by 
placing the hook on the upper grip into the eye at the lower 
end of the braided wire. The slotted bar rests on ledges 
placed on either side of the frame. The weight pan is attached 
to the lower grip and the jack 1S raised to support the pan. 
The required weights are then put in place on the pan. The 
telescopes are adjusted so that, with the crosshairs centered, 
the crosshair in the upper telescope is aligned with the upper 
scribe mark and that of the lower telescope with the lower 
Scribe mark. The test is begun by releasing the jack and 
allowing the weight to be applied to the lower grip. As the 
Specimen strains the scribe marks move avay from the cross- 
hairs. The amount of elongation can be measured by keeping 
the centered crosshair in the upper telescope aligned with 

the upper scribe mark and moving the lower telescope crosshair 
to realign it with the lower scribe mark. The amount the 
lower crosshair is moved can be read from the vernier. The 
scribe marks, although quite thin to the naked eye, are much 
wider than the hairline when viewed through the telescope. 


The accuracy of measurement must therefore be considered as 
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about +.0002" rather than within .0001" which is the 
graduation of the vernier. 

By measuring the movement of the scribe marks over a 
period of time the deformation history of the sample under 
the given load can be determined. Performing a series of 
these tests at different loadings yields a family of curves 
showing the deformation history at different stress levels. 
By measuring the slope of each curve at some given time one 
can obtain a strain rate, for that time, at each stress level 
tested. Plotting strain rate versus stress on a log-log 
plot will, in the case of a material whoSe creep properties 
conform to Norton's law, yeild a straight line from which the 


creep constants K and n can be determined. 


oa RESULTS 

Strain versus time for the 8 lb., 6 lb. and 4 lb. plates 
map lotted in Figures 1.18 and 1.19, Figure 1.20, and in 
Figures 1.21 and 1.22 respectively. The slopes of these 
plots are measured at times t=l hour and t=5 hours. These 
two times were selected as t=l hour is a point apparently 
well within the primary creep region where deformation rate 
is changing rapidly; and, by constrast, time t=5 hours has 
an appreciably smaller and nearly constant creep rate more 
characteristic of secondary creep. Creep rate versus stress 
memplocttcad on a log-loq scale in Figures 1.23 through 1.25. 


These data indicate appropriate creep constants as follows: 
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8 lb. plate n=3.7 Ke= 75.56, x 10 Gti hr. 
ae eee ior ante S hrs. 
~18 
6 lb. plate n=5.0 K = 4.10 x 10 @ t=l hr. 
ne eens a 25 hse 
~15 
4 lb. volate n=4.0 K = 8.30 x 10 @ t=] hr. 
m= 9.80 6 lao ees See 


It is interesting to note the marked differences in the 
creep characteristics of the samples cut from the three 
different sheets. All three sheets were made from lead of the 
same commercial grade and one might have expected them to have 
exhibited nearly identical creep characteristics. Such was 
not the case, however, and the differences are too great to 
be accounted for by exverimental Ganon OF data SCaccer. 
Gifkins (22) made a careful study of the effect of alloying 
on the creep of lead using lead-thallium alloy. He found 
that, in very dilute alloys, a change of composition of as 
little as .01% thallium would result in a significant change 
in arain size and creep performance. Creep characteristics, 
therefore, certainly in the case of lead alloys and probably 
also for other materials, seem to be quite sensitive to small 
changes in the composition of the material. This means that 
calculations based on assumed creep constants may be greatly 
in error unless the constants have been derived from tests on 


exactly the samc material as that being considered. 


39 


~©®@ 
- 


. 





500 
LENSIRE FESTS 
8 LB. PLATE 


x 
& / 
O) 
6 
0400 
Q 
vy 
LY 
1°) oo 
* QO 
wu J 29 
= Y a 
4c ay : 
—" 0300 A 
w Rey 
S 
ae a20P2 
\S YS Oo” < 
al 
ae US 
a 
0200 | 
OLOO 
O 2 c 4 
TIME (HRS) 
FIG. 1.18 


40 


00400 


0300 


w 
0200 


TENSILE TESTS 
8 LB. PLATE 





TIME (HRS) 
FIG. I.19 


Al 





0400 


-0300 


0200 


TENSILE SVESTS 
BUSCA E 





TIME (HRS) | 
FIG. 1.20 


42 





« ws 





se 
. é 
a 
IP 
oe TENSILE TESTS 
4 LB. PLATE 
as S\ ss 
0300 o 
20“ 
ee 
ers 
Z|= 
” 0200 . : 
8 cs aor 
e = 
@ on’ - Ape 
ot se 
52! Th o 2 6IOPS! ® a 
k 90%? " 
0100 7, 
GC 
2 3 4 5 6 c 
TIME-(HRS) 
FIG. 1.21 


43 





€ (INZIN) 


.O600 


-O500 


.0400 


.O 300 


-O200 
O 


TENSILE TEST 
4LB PLATE 
ao =!000 PSI 


S. 4 o 6 
TIME (HOURS) 


iGuilkae 


44 





é€ IN/IN/HR 


8LB. PLAT 
.O| : 


t=] HR 


n=3.7 


K=5,.56x10'4 € IN/IN/HR 


t=5 HRS. 
N23.7 
OO! : Kae 3xlOee 
.OOO! 





I 2 3 5 6 789 2 3 4 5 678910 


sTREss (73') sTRESS (Ts) 


FIG. 1.23 


Log-Log plot of Tensile Test Results, 
Seto ce lAts 


45 





€ IN/IN/HR 


6LB PLATE 
Ol 


€ IN/IN/HR 





t=5 HRS. | 
.OO| : n=5.0 y 
K=1.7x10'F £ 
OOO! 
4 3 4 6 78910 
PSI PS| 
STRESS (= | STRESS (ss 
FIG. 1.24 


Leog-Log slot of Tensile Test Results, 
Glo o lace 


46 





€ IN/IN/HR 


e 


4LB. PLATE 


Ol 
t=1 HR 
n= 4.0 
K=8.3x10!° é€ IN/IN/HR | 
t = 5 HRS. 
n=4.0 ' 
revel K =3.5x10° 
.00| 
0001 OOO! 
y 2 S 4 5 6789 2 3 Sein 6 78910 
PSI p 
sTRESS (751) sTRESS (75°) 


FIG. 1.25 


Laz-Loe ~lot of Tensile Test Results, 
foo La eS 


47 





CHAPTER II 
Theoretical Analysis 


me ASSUMPTIONS 

Consider the case of a homogeneous, incompressible, 
MeOLLOpiCc, Circular plate subjected to a time-independent, 
lateral, hydrostatic pressure. The plate 1s assumed to be 
sufficiently thin for the problem to be treated as plane 
stress with oe O. We assume the creep behavior of the 
Material may be modeled in the manner discussed in Appendix A, ~ 
which for this case gives: 


E, = ae lcei sgn 0, - z|o,| "sgn O 


m 
It 


5 BST [ia lhe Sgn O, t\o,|" sgn 01] 2.2 


where sgn is the signum function defined as 


ey +1 x>0 
=e -~1 x<0 


in this form 2.1 and 2.2 can be used for negative as well as 


positive stresses. 


B. BASIC EXPRESSIONS FOR MOMENTS AND IN-PLANE FORCES 


Sony MiG elerancde2. 2. TOr 0. and Op in terms of the strain 


rates we find: 


_ pen IWAN TS eas Saar a We 

oe = (Sy) ewer sgn (e +2€,) 2.3 
Meelis) L/n aan 

an (3%) oneal | sgn (€,+27€_) 2.4 


The extensional strains along the middle surface of the 
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mieate will be denoted as e.. and Cg et omemeratuce of the 


neutral surface as Se and Kg: Adopting the Kirchoff assumption 
that sections which are plane before deformation remain plane 
during deformation, the total strains, E and €,, at any point 
along the thickness of the plate may be written as: 


€ =e + Zk | 2.5 
r r r 


Eg = e amok ZO 
where z is the distance from the middle surface. Notation and 


Sign conventions are as shown in Figure 2.1. Taking derivatives 


of 2.5 and 2.6 with respect to time gives the strain rates: 


See= ©. + ZK Deco 


Meet itutinag into 2.3 and 2.4: 
7m 


ee 4 ya * ] « e ics 
oS are |e. GAS Z(K_+2KQ) | 2.8 
« le e le 
sgnle,+7e, tZ(K +3KQ) | 
= 4 ,1/n @ le e le Jiy/Aig) 
Ba (35) Qe +P lel ae Z(K,+2k) | Ze 


: 1s . le 
sgn [é,+76, +2 (Kytzk_) | 


Peon this point conventional notation for the signum function 
Will become rather onerous and we shall define the symbol 


Bane Pp 
J£(x,, Xx el 


pe 
as ‘meaning the absolute value of the function f is raised to 


the power p, the result is then assigned the sign which 


fenctkion f(x, x 


5! ae assumes when evaluated with the 
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appropriate arguments (Xx), Xo eaeae ) for the case in 


question. 


Considering a homogeneous thin plate we may compute 


in-plane equivalent forces per unit length, No and Nas as: 
h 
N= 2EOndZ 2.10 
r 
fh 
bo 2 oe 2.11 
e 
Equivalent moments per unit length are: 
h 
_ 2 
M = | Oa 27 ieee 
r E 
2/8 
2 
A 
M, = 2 0 9zdz 2.13 
-) 


By substituting (2.8) and (2.9) into the above equations, 
fet0) through (2.13), and integrating over z, we may obtain 
expressions for the plate force and moment resultants, Nous 
Nos Me and Mg: Some care, however, must be exercised in 
performing this integration. The signs of various terms in 
the result depend upon the relative magnitude and directional 
sense of the extensional strain rate and the curvature rates 


in the plate. It can be shown, however, that upon integration 


all cases may be represented by expressions of the form: 


SAL 
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re ctoeey I 
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0a S) 
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7 sta yee 
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oe Te i n 


* eee 
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ce. DEFORMATION PROFILE AND STRAIN RATES 
For a circular plate undergoing finite deformation the 
curvatures and extensional strains can be written: 


oe es! 
. = —WwW Kp = 2.18 


E ae 
= i A a 


where w and u are the transverse and radial deflections 
respectively with directional sense as shown in Figure 2.1. 


Taking derivatives with respect to time we obtain the rates: 


age 


De 

II 
S 
+ 
=e 
=> 
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il 

HIE 

Nh 


To proceed further, the transverse and radial deformation 
rate profiles for a plate with a given set of edge conditions 
must be determined. For the Biplive IAaeee elastic deflections 
of a clamped circular plate Timoshenko and Woinowsky-kKreiger (23) 


employed the elastic bending analysis profile: 
v= Ww (y-Eoy2 De 
O a . 


where We 1s the deflection at the center and a is the outer 
radius. In the experiments described in Chapter I it was 


moved, for the = ratios examined, that the measured deformation 


profile could be represented by equation (2.22), especially at 
W 
lower values of ma Assuming, therefore, that the vertical 


deformation profile may be represented by (2.22)the curvatures 


and curvature rates may be written: 


ANS Saee ae xr? 
o 7 are Kg = —az ing?) pee 
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: : ; 
a e ee Orne As 
ees -7 (l--—7) Kg a lee) 2.24 
For the deflection rate in the radial direction, t, we 
adopt an expression suggested by Timoshenko and Woinowsky- 
Kreiger (24) for large deflections of clamped circular 
elastic plates which, adapted for time dependent problems, is: 
2 


u = r(a-r) (c, ae Oe) Ai C3r to ae 


n 
5 1 aay Te Zee 


oe 
This expression satisfies the end conditions u=0 at r=0 and 
r=a and provides for a wide variety of possible profiles 
through variation of the Cc. values. The constants, Cis will 
be determined by requiring minimum internal energy dissipation 


in the plate. Taking the first two terms of the series we 


have: 
7 ; 2 ae 3 
iw = c,ar c,r + cjar Cor 226 
from which: 
&' = (a-2r)c, + (2ar-3r°)c, 2.27 
a eave wee (ar=r-)c | 2.28 
r 1 2 ° 


With these expressions the extensional strain rates, (2.21), 


become: 
16w we 2 
— _ eee, OO ee. 2 
e.. = (a 2r)c) + (2ar-3r )c, +> Saat ee 2.29 
>: = feo. “eee 2.30 
¢) il Dg 


Examining the expressions for é. and e 5 we are led by 
dimensional considerations to suspect that Cy and C5 are of 


the form: 
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SY ee ie ere, ee ore) 

a la ey ae 2.31 
where Cy and Co are pure numbers. In order for each term in 
the sum to have the time dimension each must contain Wo If 


ag and ep are to be zero at any point on the plate when 


a 0, then each term should contain Wo To be consistent 


with the last term in (2.29) the diameter, a, is selected as 
a length factor. If one considers extensional energy 
dissipation only and solves analytically the n=1 case for Cy 


end Cc it is seen that they are indeed of the above form. 


of 
ae Bees , 
Defining x = 3, We may now write the curvature rates and 
extensional strain rates as: 


: = 2) Dele 
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A 
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: 
a 
1 
ee) 
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Ae 
il 
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a 
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Ro 
— 
NO) 
Ww 
WW 


é. = Se (lees ee ee elon =x) =] 2M 


WoW : 
4 an [(1-x)c} + (x-x*)c,] 2.35 


(De 
II 


We now have expressions for the curvature rates and exten- 
Sional strain rates in terms of geometrical parameters for a 
Clamped circular plate. These expressions can now be sub- 


stituted into 2.14 through 2.17, the force and moment resultants. 


W 
For the higher values of — the experiments showed that the 
deformation profile tends to the parabolic membrane form: 
2 
> e4 
a: We ‘ea aoe 2.36 
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The problem may be solved using this profile also and, as 
will be seen, the solution for the center deflection rate, 
as tends to converge to the same value for the two profiles 


Ww 
as = increases. The same procedure can be done with other 
assumed profiles as seen in Appendix D for the simply 


Supported case. 


Pee ENERGY EXPRESSIONS 

The total internal energy dissipation in the plate during 
creep deformation may be expressed as: 

’ 271 7a 

V= | | (MK Lt+M,K tN € +N ,6,) rdrdé 2238 

0 0 

where transverse shear energy dissipation is neglected. The 
integration, strictly speaking, extends over the deformed 
surface area of a plate. However, if attention is confined to 
moderate transverse deflections and small strains then the 
deformed surface can be taken to equal the original surface 
area. This procedure is usually adopted for deriving a 
consistent set of equilibrium equations and strain relations 
when using the principle of virtual work. (25) Using the 
pre, tously derived relations, equations 2.14 through 2.17 and 


2.32 through 2.35, we may now write for the clamped case: 
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eee 2, | / mien 2+1/n uw ky Aor Gay ilar dl ye 
Beee “SK? aot 2 (3) ap (IE+D| 
W W 
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With these expressions 2.37 can be written as: 


W 
= Dy a Te 2+1/n o,1l+l/n 
SS iia) are so (7) 
1 
i (is + ie + fA ” fy )xdx RNY, 
8 a6 8 


where the quantities ie are call functions Of x, Cys Cc 
b 


ma and n. These symbols, f. , represent the terms containing 
b 


xX in 2.38 through 2.41 which must be integrated over the radius. 


The result of this integration over x will be a function of 
W 


—, and n which will be denoted I. The equation for 


no? h 


ie” 


V may now be written as: 


W 


CHG) 243 


4 ,l/n n y2ti/n 
il 


Wo 1+1/n 
(=3) 


The quantities Cy and C5 aS previously discussed determine 


the radial deformation rate profile. They must be chosen such 
that the rate of energy dissipation in the plate is minimized. 


Clearly, minimizing V with respect to Cy and Cy is equivalent 


to minimizing I with respect to Cy and Co- The minimization 


of I must be done for the current value of — and n. The 


correct values of Cy and C4 for one set of values for the 


parameters ae and n may, in general, be different from those 


for some other set. One way of determining the appropriate 


values of Cy and C5 might be to substitute values for — and n 


into an expression for I and take derivatives with respect to 


C7 and ci, set them equal to zero, and solve for the values 


Oe 
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resulting in a minimum I. Unfortunately, the expression 


1 
i (fi + se + i + i ) xdx 
r 8 r 9 


is such that the integral cannot be evaluated analytically 

so there 1s no analytical expression for I which can be 
differentiated with respect to on and Co. We must, therefore, 
evaluate I numerically and search for those values of Cy and 
C5 which result in the miminum value of I for a given case 
(wo/h, n). This was done, for the cases considered, with the 
computer programs whose listings are given in Appendix H. 
These programs were used with an IBM 1130 Computing System 

at Massachusetts Institute of Technology. The computer 
output yielded a series of smooth curves for regular variations 
of Cy and ci. These curves showed clearly identifiable 


Z 


minima with respect to Cy and Cc. for the function I. The 


procedure is illustrated for a typical case in Appendix E. 


Oe SOLUTION FOR DEFORMATION RATE 

If a structure is assumed to deform in a manner prescribed 
by a kinematically admissible virtual velocity field, (i-.e., 
continuous and satisfies the velocity boundary conditions), 
then the principle of virtual velocities states that the total 
Virtual internal energy dissipation is equal to the total 


virtual external energy dissipation. 
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Ww 
Having minimized the integral I for some case (—, n) we 


find the deformation rate, Wor by using the principle of 


virtual velocities and considering a small variation of the 


Geformation rate about its equlibrium value. This gives: 
av - = s rer 
aw. Swe = 27 J qé6w rdr 2.44 


where gq is a time-independent, transverse pressure. For the 


case of the clamped plate we have: 


2 2 
2n f“qéw rdr = 2nqéw_ f*(1-=7)?rar = “22 ow 2.45 
0 O 96 a 3 Oo 
and: 
Aum, . 1/n Ome ler re 2+ 1 /ni as ~. 
aw °“o 7 [21a (3%) 2) h T]ow, WARNS 


The expression 2.46 is the derivative of V for some specific 
W 


value of == and the derivative is properly taken considering 


W 
O ; 
A as a Seiisednteeewe are seeking the value of Wo for a given 
W 


instant when a= has some fixed value. An instant later Wo 
will, of course, have some other value and at that time Wo 
will also have some other value, but at a given instant wi/h 


may be considered fixed. Equating 2.45 and 2.46 we obtain 


the result: 
qzt2/n 


ck n 
Wo 7 A/mm 241 /n Beto 
(35) h 61 
or 
Yo o= [ (Sy 241 /n (3Kay 1/n cy)” > 48 
O h a TET 
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If we define the quantity ©? as: 


W 


(——) 

Kagn Wag 2.49 
a, 2n+l i 

3/4 (>) 


mas Can be written in the non-dimensional form: 


W 
O _ nee figek n 
(Kaqi? = 3/4 (F) oO. 2.50 


The parameter ? 1s a direct function of I, the form of this 
function depending upon the assumed velocity profile. For the 


clamped case: 
ee: 
d = (FF) Bao 


Tt 1S convenient to define ¢ as it permits a direct correlation 
of results among the different velocity profiles considered. 


When the quantity I has been minimized with respect to Cy and 


W 


Cor @ is a function of n and = Values of 6, for the clamped 
W 


case, are plotted versus ae for various values of n in 
Figure 2.2. Use of these values with any of the equations 2.47 
through 2.50 will give the center deflection rate we as a 
function of the applied pressure and the plate's material and 
geometrical properties. Once we has been determined, the 
deflection rate at any other point on the plate can be found 


from the assumed velocity profile, which in this case is given 


iby equation 2.22. 


61 








U 


[o* 


a 





ec Old 


MT) anal) 


$3903 Q3dWV19 
3LV Id YV TNOYIO 











62 


It was noted in the plate experiments that for the 
W 
higher values of oa the profile began to tend toward a 


parabolic form (or membrane profile) described by the 
equation 2.36. The entire analysis just described can be 
repeated using the profile 2.36 vice 2.22 to obtain values 

of ¢ for the parabolic configuration. When this is done it 
is found, as shown in Figure 2.3, that the difference between 


W 


the two solutions is small for —— = 1.5 and has virtually 
disappeared when x = 3.0. It seems, therefore, that when 
the plate deforms to a value of = large enough for there to 
be a significant deviation from the profile of equation 2.22 
it will have deformed enough so that the solution for We is 
relatively insensitive to the assumed profile. 


The computer program listed in Appendix H is written for 


an assumed velocity profile of the form: 


; : ne Ta? cs 

w= ent a—7) 4 a5 (1-=z) J 252 
where a, + a5 = Pee ne analysis ueilizing this profile is 
outlined in Appendix C. If we assign the values a, = 1 and 
a5 = 0, the analysis with 2.52 is the same as for the clamped 
profile,2.22. If we assign values a, = OQ and a. = 1 the 
analysis is for the parabolic membrane profile, 2.36. Some 
experimentation with fractional values for Qa and a. was 
attempted, but these vroduced no more useful results. The best 


aqreement with experimental data came with Oy and a5 at their 
extreme values as above. The final result using this profile 


this 
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Coy™varison of clamsed and rembrane solutions 
tone — 4,0 
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js a 
(2 27 ae Cee ae 2.53 
n h iZ 
Kaq & 
SO in this case 
3-0, 


where I. is the appropriate integral for the combined profile. 


ere DIMPLY SUPPORTED CASE 

The simply supported circular plate may be treated in 
the same manner as above if the velocity profile appropriate 
to it is known. In the elastic bending analysis of this case 
Timoshenko and Woinowsky-Krieger [26] determined a deflection 


profile of the form: 


Pa 2 
om 21g es 
Beech eae 051 a7 oo 


which for v = 1/2 gives the velocity profile: 


me ae Seg 
Ww = w, (1-22) (1-7, ey ° 2.56 


There is no experimental data available in order to verify 
this profile as in the clamped case, but at least a reasonable 
estimate of the creep behavior of the simply supported case 
can be expected from an analysis with the profile of 2.56. 
This analysis is outlined in Appendix D and yields a final 


result of: 
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so that in this case: 


= 5 
¢ = (S577) 758 
S 
where Te 1s the integral appropriate to the simple support 
profile. Curves of ¢? versus — for various values of n are 
given in Figure 2.4. In Figure 2.5 the results of the 


clamped and simply supported cases are compared for the 


case n= 4. 
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CHAPTER Lit 


Comparison of Theoretical and Experimental Results and Results 


of Previous Investigations 


A comparison of the experimental and theoretical results 


is shown in Figures 3.1, 3.2, and 3.3. Here the non-dimensional 


W 
quantity (-2~) /(2)?"*? is plotted versus the non-dimensionalized 
Kaq W 
center deflection, — for the creep parameter n appropriate 


to each of the three types of plate tested. In these plots 
the appropriate analytical results for comparison with the 
experiments are, of course, those for the clamped case. Also 
shown are plots of the bending analysis results of Venkatraman 
and Hodge [1l], and of Odqvist [12]; and the membrane analyses 
eepedqvist {12](16], and of Onat and Yu ksel [17]. The 
Mathematical details of xeducing these various analyses toa 
form appropriate to this plot are given in Appendix B. The 
bending analysis solutions are horizontal lines Since that 
type of analysis does not consider geometry changes. The 
membrane analyses consider geometry changes but not bending 
moments within the plate. Experimental points are given at 
each pressure tested for times of 1 hour and 5 hours. These 
points are found by obtaining the deflection rate at the 
indicated times from the plate test experimental data and 
substituting this value along with the appropriate creep 
constants found in the tensile tests and the pressure and 
geometrical constants into the form (2 yf (2) ont, The 


Kaq™ 
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FIG. 3.5 


Comparison of experiment and theory, 
Apmis plate 
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value thus obtained is then plotted versus the appropriate 
WW 
value of — as determined from the plate test data. The 


h 
times t=l hour and t=5 hours were chosen as representing 
cases of rather rapid and relatively slower deformation 
respectively. 

The marked influence of finite deformation on the creep 
behavior of a clamped circular plate is clearly indicated by 
both the theoretical prediction of equation 2.50 and the ex- 
perimental results, which are seen to be in substantial 
agreement. 

A similar plot, Figure 3.4, shows the relationship of 
the simply supported case analvsis to the results of the 
Simple support bending analysis of Venkatraman and Hodge [11] 
and of Odqvist [12] and the membrane solutions of Odqvist [12] 
{l6] and of Onat and Yu ksel [17], for the case n = 4. No 
experimental data is available for this type of end condition, 


but the analytical results also indicate a marked influence 


from finite deformation. 
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Summary and Conclusions 


Experimental and theoretical results have been presented 
which show that the effect of finite deformation upon the creep 
behavior of circular plates is of primary importance. It is 
necessary to consider this effect even when the maximum 
deflection of the plate might be well below the plate thickness. 

A method of modeling the constitutive behavior of a 
material undergoing creep deformation under a multiaxial state 
of stress has been proposed. A technique for solving multi- 
axial creep problems by using an energy method and the 
suggested constitutive equations yields a result in the form 
of an equation and graphs which eta be quite easily applied. 
Solutions were obtained for a circular plate under a time 
independent uniform lateral load with clamped and simply 
Supported edge conditions. A oe ee of theoretical and 
experimental results for the clamped case showed substantial 
agreement; indicating that the analytical technique described 
may be expected to provide quite reasonable estimates of plate 
deformation rate when the creep behavior of the material may 
be adequately represented by Norton's law, and the appropriate 
values. of the material creep parameters are known. A comparison 
of these results with those of several previous investigators 
has also been presented. Lead was selected for the experimental 
material in these tests but the analytical results should 
remain valid for any material whose behavior is adequately 


represented by the constitutive equations. 


tis, 





Specific consideration was given only to circular plates; 
it can be expected, however, that the conclusion as to the 
importance of finite deformation will apply as well to plates 
of other shapes, initially straight beams, and to initially 


Straight cylindrical shells. 
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APPENDIX A 


Coistititrve Equations for Multiaxial Creep 


Creep data is not plentiful, and that which is available 
is almost always the result of uniaxial tests. It has been 
found that much of this experimental creep data can be repre- 


sented by an expression of the form: 


e = Ko” (A.1) 
where K and n are material parameters which are, in general, 
functions of time and temperature. This expression is 
commonly known as Norton's law. [3] 

In generalizing Norton's law to a form suitable for 
multidimensional problems the most common approach has been 
to use the concept of equivalent stress and equivalent strain 
rate for determining the proportionality between strain rate 


and deviatoric stress. Mendelson [27] uses a relation of the 


EyDe: 


Me 


Pas (A.2) 
where creep equivalent strain rate, EQ! may be related to the 
equivalent stress, Oo! by Norton's law. 
Pere : 

U, = 35, (A.3) 
and zs is the strain rate resulting from a stress of this 


magnitude. Odqvist [6] uses the expression: 


Og Shee 
Syed 1) 
"aay 2 eo OK (A.4) 


where E is some standard value of creep rate for which O=0.- 


dey 


These expressions seem reasonable and satisfy the mathematical 
requirements for multiaxial creep constitutive relations, but 
are often rather intractible for analytical work since one 
must always compute O, in order to know the proportionality 
between the local stress and creev rate. 

In this analysis the constitutive behavior of a material 
undergoing creep deformation in a multiaxial state of stress 


Will be modeled in the following way: 


a es, ee RAs n 
eT Klo. v us Vv oy] (A.5) 


where the subscripts 1, j, and k represent principal directions 
of stress and strain rate. It can be shown that for conser- 
vation of volume we must have v = 1/2. In real materials such 
things as voids and inclusions may give rise to situations in 
which volume is not conserved. If this effect is important 
then in such cases some other value for v must be assigned as 
appropriate to the specific material under consideration. 

It is also desirable that the constitutive relation (A.5) 
hold for both positive and negative o for all values of n. 
The significance of positive and negative signs here is merely 
One of direction, tensile or compressive, and this directional 
sense must be retained for all values of n. This may be 


accomplished through the use of the signum function defined as: 


{ +1 x >0 (0.6) 


eee ta x <0 


An expression of the form |x|" sgn x is used to denote that the 
absolute value of x is raised to the n power and the result is 


assigned the sign of x. 
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For homogeneous isotropic materials, therefore, the 


peoeosced Constitutive relation takes the form: 


e.=K[|o, |” sgn o,-2]o,|" sgn o5-2] 0, |" sgn 0, ] CRIS 


These equations are valid for a material which has the 
following properties: 
(1) Homogeneous, incompressible, isotropic material. 
(2) Coaxial stress and strain rate tensors. 
(3) Zero deformation rate under hydrostatic pressure. 
Pymecercinates correspond to principal directions of 
stress and strain. 
(5) Creep behavior corresponds to Norton's law in 
uniaxial stress. 
Lmanoubhemocmotca that this model predicts an effect from 
a hydrostatic pressure superimposed on a deviatoric stress. 
The validity of this may be questionable and was not really 
tested in the cases studied in this investigation. ~ For 
applications where larger hydrostatic pressures are involved it 
may be better to use only the deviatoric components of the 
meGesSSes. 


PEE memedseeOr ad thin Circular plate, choosing cylindrical 


coordinates and assuming ote aCe O, equations A.7 become: 
_ = KE ]o.|- sgn o, =a |e sgn Oy] 
E a = Kop |” sgn 0, phe sgn J] | (A.8) 
E. = Sey oo is sgn o,. -t|o,]” sgn o,]. 
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APPENDIX B 


iieernretacton of Results of Previous Investigators 


Introduction 

The following theoretical treatments all apply to circular 
plates loaded with a time independent pressure uniformly 
distributed over the entire plate area. The infinitesimal 
solutions in Sections I and II retain only bending moments 
while solutions in III and IV consider geometry changes but 


retain only membrane forces. 
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I. Venkatraman and Hodge Solution [11] 


a. Clamped Edges 


The maximum deflection rate, at the plate center, is 


expressed in the form 


w= wt eet a 
Oo 


where G iS an arbitrary function of time and 





u= 2A 


2n+] 2 


The term \ is obtained from the assumed form of Norton's law 


given as 


e = (>)". 


n (By 2+1/n_ 


(Be. 


CB 


When compared with the form given by equation A.1 which was: 


a) Ko 
we have 
= . l/n 


Reeereleucing B.1.2, B.f.3, and B.1I.5 into B.I.1: 


| ere Gn Zigarle 
Wo W G(y) 


OF 


a, 2nt+l 
(2) 


Kaq 


nN 


92+1/n 


[ 


24n 


(2ntl) jn 


(Bi. 


(B.1. 


(Bi. 


(B.1. 


which interprets B.I.1 in terms of parameters introduced in 


Chavter III. 
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aA) 


2) 


3) 


4) 


5), 


6) 


7) 





For comparison with the results in this paper we take 


G=l. The term w* is plotted as a function of n in reference 


ect). 


b. Simple Support 


For simply supported edges the deflection rate at the 


plate center is of the form: 


Q 
We = he (ae (a CE (B.1I.8) 


oa) 
n 2+1/n 
Kaq ae 2 (2n+l),n 


a, 2n+l 
(2) 


where as before w* is plotted in reference [1l]. 
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mie Odqvist Bending Solution [12] 


This result is of the form 





We =al™ v(o) (B.II.1) 
where ee 
_ Pageell Sele qa 
it = we ho. (Bakr 32) 


and V(o) is plotted as a function of n in reference [12] for 
the clamped and simply supported cases. The term oF 1s from 


the assumed constitutive relation: 


oO Shy 
>, SS ee ee 
€1j 20 O CBr 3) 
Cc Cc 
which for uniaxial stress reduces to 
° 1 n 
e€ = Oo n 0 (B.II.4) 
Cc 
and when compared with our notation: 
€ = Ko” (B.I1.5) 
gives: 
meee. Ln 
a - (z) (Be 11.26) 


Substituting these expressions into B.II.1 and rearranging we 


have: 
Wo 
a ASG SSE 9) wea 
2 (a) (3) 2 V(o) (Edi?) 
aventL a 


oS 





III. Odqvist Membrane Solution [16] 
In this analysis the transverse deflection at the center 


of the plate is expressed in the form: 





Zo" h Wo) 
ee (pes ii) 
O q 
where 
2 
© 24 “ 
W(o) = Te (1+gc+ (VR eet 2) C2}, erie) 
ee OE 
- Sra, (Bevin 3) 
Cc 
= -C, B. ; 
8 (n+3) 50? 
1S 
a [=p + 8(n+2) (n+3) oe) Bes (B. III. 5) 
s 0 
O 1 


h is half the plate thickness, and C is tabulated as a function 
of n in the reference. In this form equation B.III.1l1 gives the 
center deflection as a function of time. For the deflection 
rate we take the derivative with respect to time: 


; 20 hi . : 
v= z, TI) (B.III.6) 





Potferentiating B.III.2 with respect to time: 


(An oy hn JE 


0 e 
; es s 
W{o) = r 2 {1+ 4c + 18 tats) (B.III.7) 
from B.III.4: 
: eto 
g = -c Biot3)s (B. III.8) 
Z 
0 
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Mist lewtang &.111,5 and B.-1II.8 into B.III.7: 


: 2 2 
w(o) = 2C(n+3) [1+4c+ f4n7_+27n + 9)C') 





18 (n+3) 
n+l 
[ney + 8(n+2) (n3) SE ne? (B.IIT.9) 
S. e 
meester eutLing B.J1I.5 into B.II1I.2: 
2 
Pp 1 
eee et ct j~—> 
wlo) = n [ Bo te) nt 3) ~|n+2 
2 e) 
O l 
2 
(ile Mow eens ocala DD ne C2] (Beckie) 
polving for t: 
2 
a 4 wlio) ,nt2 1 
© * ectnray mea) |B) na! eee? 
a O 
= 1 (4n? +27n + 9) 2 
where B= 1+ 3C + Sa ees Cc : (Be iit 212) 
puostituting for t in B.III.9: 
DN ee 
So) = 26 Bor” _{nt3) , 
+ + 
Uo (0) )™*+ 
0 2n+2 
; 1 
> 
(wio)) ttt (B.ITII.13) 
where 
5 . 2c B™T* (n+3) 
n+1 : Beeps 14) 


4 
Note that 6 here is that of reference [16] and bears no 


relation to the quantity $ as used in Chapter II. 


85 





Odqvist's constitutive relation was: 


Oe 


= 3 (277 4) (Beit. 15) 


aa (=) 9 oD (B.III.16) 
e 
iyi ng : 
eae 
as Le (Beli. 17) 


eiieattewioamih. |, Bolii.3, B.LTIL.6, B.III.14, and B.III.17 


into B.III.13 and rearranging: 


Wo 
. Kaq™ | : 2h,ntl | 
~_  e ae = 96 i) (B.ITII.18) 
(a_y 2ntl 2 
2h 


Odgqvist uses the symbol h to denote half of the plate thickness 
rather than the full thickness. If we call the full thickness 


h, we have a form analagous to that of the other results: 


Wo 

< 2 1 
ELS Ee 

(a) envi a) Ww =e (B.II1I.19) 

h Sime 
where 

* + + +r 
pee 2c{1+4c + “a Cele = (n+3) 
antl 


and C is tabulated as a function of n in reference fl6j. 
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IV. Membrane Solution for a Circular Plate, Onat and Yuksel [17] 
Onat and Yuksel assume that an initially flat circular 

plate deforms into a spherical surface upon application of a 

lateral hydrostatic pressure. Consideration of this geometry 

shows that the middle surface strain rate tensors have prin- 

cipal directions which are directed along the sphere's 

parallel, meridian, and the normal to the surface. The first 


two strain rates are positive and the third is given by 


= a = eee CO — ' a 
a ae atan 5 (Bo 2) 


For a spherical membrane: 


D 
2 ee en OO fe Pe, 
Tie - o oh dh. Sina (itcosa) (B.IV.2) 


Creep deformations are governed by the relation: 


a tan 5 = Eo) (Bery. 3) 


where for Norton's law 


f(o) = ¢€ (—) ee OS Be, (B.IV.4) 


The symbol p in B.IV.2 above represents the uniform pressure 
per unit area as used in [17], dD, is the diameter of the built- 


in support, and he is the original membrane thickness. 


she 





Pecaine Bays 


From the geometry of the above figure we see 


w_ = R(l - cos a) 
e) 
6 
7 = Rsina 
Sin a we 
Egquating R mie cose oo 
a ; Cae aCOse 
By trigonometric identity: tan > = Saga 
O avs 
therefore, tan 5 soe 


Differentiating with respect to time 


d a O Oo 2W 
=~ (tan=) = = sec’®’= = —— 
he 2 2 2 D 
O 
4w 
or a =—2 
D sec?& 
Oo Z 
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(Backy 


(Baby. 


(Bev. 


(Beay. 


(Be Ly. 


5) 


6) 


7), 


8 ) 


9) 





eae Nom bevel eb. V2.4 ,05.8EV.8, and B.IV.9 into B.IV.3 





Aw 2wW Dae 
(——-—) ( =“) = K (—?—)"” a 
D 20 D Zia Sina (1+cosa) 
sec*> O O 
O 2 
and rearranging 
Wo 
(——+——) 
Kap 
——————— (a sec?=] 
W Sina (l+cosa) 2 
a,ntl O 
(>) (2) 
O O (B.IV.10) 
wiiere a= Po : 
2 


The non-dimensional parameter resulting from the energy 


analysis 1s 








Me 
( =) 
Mee ee, dn zt 
(a) entl - (ayn Ww, ener (ltcosa)’ secs] (B.Iv.11) 
oS hy (=) 
O 
Wo y 
Given a and Te angle a can be found from B.IV.7 and values of 
2 O 
W 
O 
Kap we e 
a. 2ntl1 computed as functions of a and n° 
(=) 7 ‘ 
O 
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APPENDIX C 


Mathematical Detaiis for Analysis with Combined Bending and 
Meneaanewe morte 1Or the Clamped Circular Plate 


Experiments indicate that a clamped circular plate 
undergoing creep deformation from a transverse hydrostatic 


pressure assumes a velocity profile described by: 


2 
° = e ae 4 
oo W 

when i 1s relatively small. At larger values of — the 


orofile tends to the membrane form given by: 


: : 32 
= ek 2) (C72) 


These two forms can be combined by writing the velocity 


profile as: 


oo 1D ag 

Ww = Wola, (1 52) at a, (1 a (C3) 
where 

Cy + a. = ] (C24) 
and: 

0 1S es a 

0 < a. s 1 (C25) 


Paudatton C.3 sredtices £6 the small deflection profile, C.l1, for 


a, = aoa ees one €o the membrane promi le, C.2, for a, = oe 


a, = 1. Using this profile and proceeding as in Chapter II: 
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° e e 2 

oe = -w' = wi loy oa (1-277) + a, a (C.6) 
2 eee) + a, 25] (ena) 
a “y az" "az a ; 

and ; 

: ; we r2 

Sas 2K, = se [[20) (Steve) 2 2 (E28) 
; : Wo 2 

Ky 4 2K = aes [2a, (3-527) + 305] A (C.9) 


In these equations and those to follow we may, of course, 
eliminate either a, Or a, by using C.4. This was not done 
here, however, in order to keep the terms from the C.1 profile 
and those from the C.2 profile directly visible and separated 
in each equation. 

For the radial deflection rate, u, we use the form: 


u = r(a-r) (c,+c,r) | (CC. 0) 


The extensional deflection rates are: 


ee a - | 
e- u' + w'w eg = (Gaae0) 
Using as before; : 
r | lc Hays 
= = oc! =- oe! 
Mee 2 1 gy 52 = %2 ar 
: WoW 
a ay ' a Zee ART 20 UE | gate Sea, ee 
e.. = i ea 2x) Cc) + (2x-3x )c4 + 16x a7 (1 i) 
2 _o2 2.2 
+ 16x (Ly Ay (1 x“) + 4x a5 | (C22) 
tg, 
== = ' Ee i ae 
Eg ee { (1 x) Cy + (xX-x )cs] (C.13) 


ot 





These reduce to the avpropriate expressions for the small 
deflection profile when a, = ioe. a. = 0, and to the appropriate 


membrane profile expressions for a, = 0, and a, = i 


W WwW 
: vi O O re i Sei ; ny a a ey 
e. + WE san 3x) Cc} + (3x ZX) Ch + 16x a5 (1 x*) 


2 eine 2.2 
+ 16x aja, (1 x°) + 4x ab J (C214) 
e e WoW, 
= Bt dee 5 ae Z 6 
Sig + Ze. = apse [ (+ 2x) Cc) =P (2x 3x Co 


+ 8x7a2 


ee 22 2 Ee 2.2 
ae x°)* + 8x aja. (1 x“) + 2x a5] (c.15) 


The general expression for moment per unit length derived 


earlier 1S: 


Peevey teen, 1 eh) bhjiti/n _ n Di) 2a 7 
eer eer 5 62! 2tD | Wei vole 
¥ 
ig). Jodo dense n penn 2ti yn 
+ gla-z (ntiyp! 22 ls ° Oe 


and for extensional forces is: 


N= (=) n dL bh,1+1/n 7 _bh thepil As) 
3K Scere lene la =| (Cai) 
" : ae ieee 
Where for te and NU: a er + Zp 
b= + #£ 
1 a 9 
and for My and No : fe é, i Le 
b= ££. + #£ 
8 
Defining: 
= co me 
ae = 2a, (3 7x°) + 30.5 
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W W 
aS 


< nee » O Ba arene ' 
(3 2x) Cc} = + (3x 5X }cs - 


oe 
1 


W W 
Pay 2 Ww 2 a5 SERS) ee 
+ 16x a; (1 Kenia  Oomct 1 6X A, a5 (1 xs) A + 4x 5 
= m 2 
D. = 2a, (3 5k =) y + 30.5 
Wo ve Wo 
~ a5 » O _ 2,1 O ee ee oO 
BE. = (3 ZS C y (2x 3x Jes + 8x ay (1 ee 
ie ve 
2 Memos © revs 6) 
+ 8x aa, (1 32.<5) h + 2x a5 pe 
we may write: 
W 
ee eee t/t) 2+1/n Ogletree, 2 1 1+l1/n 
MM. ‘co (32) n+1 h eee 2A |B. i A 
2 n 2th iE Lai 
fonemmacediss  f8clg 9 > ga_! Be ~ 2A, 
Cc Cc 
n 7 Pin AD fa 8 ae, 
“(antl) Az | Bo 2: ] [4a, (1 3X )+20,] (cC.18) 
wW 
2 ere eas Zen 7 eOwl tym 1 Tepe A 
rem (3%) n+l : ae A. [|B 2 $A | 
W W 
. _ lenge ees ae eo. 7 Nei ee) 
[Bere 2a. | J ele2 ete t (2x > 3x") C> = 
Wo Wo 
22 ee ND oe =sr 2 abd 
+ 16x a4 525) h + l6a,a, (1 aan a 
Wo 
Oo ie) 
+ 4x 5 5 ] (C.19) 
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W 
eet 2+1/n Owen Zz 1 1+1/n 
haa (3x) _ ntl 2 (7) lap |B ss Da! 
7 n 1 Mande 1 _ Ty n 
(2n+1)D2 | ee Bell g op | EG 2 DI 
Cc Cc 
n ee Bee Lyi Bato 
+ TintlyDe 5 = wale ] [405 {2-x*) + 20,5] (C- 20) 
Oey an 2+h/n ,10,1+1/n 1 1, jl+i/n 
Noeg = (Sy) sar (7) Do (JE, + 2DQ| 
-le - 3p [2*1/7) . pexyot 2 + (xex2)02 -2). (C222) 
8: ac ei Oe ige. in ee 


The total energy dissipation rate is 


obtained by integrating 


these exvressions over the deformed surface area of the plate. 


This gives: ; 
1/fn n 2+1/n Wo 
Ss ey 


n+l] n 


4 
2 ee 
Dia (3%) 


V 


where I, is the integral of the bracketed 
: W 
over x for a particular case (—, n), and 


1+1/n t 


e (C222) 


expressions above 


where the integral 


has been minimized with respect to Cy and cy for the value of 
ay and a. under consideration. 
Taking a small variation of w_ about its equlibrium 
W 
position for some particular value of a! 
av st - a . 
dw OW _ 27 é g Sw dr CC 1250) 
On 
e l/n wee meen tl /n oT 
We [27a (3%) h ae T.] OW = 
2nqdw fae el eae + eid (Ge24) 
Ue Oo Oo oo az 2 at : ° 
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q2t2/n - 
Poy Rape i ep J (oee2), 
(=) h ——_ [ 
3K 3-a Cc 
1 
or 
Bee 
. _ awe Ha oka, 1/7 I n 
oe [ (FP) aac ) a ae ' (C.26) 
and in non-dimensional form: 
W 3-a 
O a) Ci otal eS Tt 
Kaq Cc 
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APPENDIX D 


Mathematical Details for Analysis of a Simply Supported 


Circular Plate 


In calculating the elastic deflection of a simply 
supported circular plate subjected to a uniform lateral 
pressure Timoshenko and Woinowsky-Krieger [26] find that the 


deformation profile assumes the following form: 





7 1 l1+vy r? 
Ww = Me (1 aoe (1 een a) CDi) 
where Wo fiemeilcamGeiceh aderlection. if we make vy = 1/2 and 


take derivatives with respect to time, we get the deformation 


rate profile: 


2 2 
3, ule eG oe 
ie? Gil ma Ti Ate (D. 2) 
One may now compute: 
W 2 
. nue Qe (2 SEG: 
eee at ii il az’ SE ey 
ae W 2 
_ 6 2 ie (2a ees 
eames s af an Tl az?’ a 
and 
Sly ye 7b) ee 
Re ty = az lip at az! — 
W 2 
be = 9 (42 30 x 
Sp iS. = => Pees Bee) 
Using, aS in Chapter II, the exvression: 
i = r(a-r) (c, ac) (D.7) 


for radial deflection rate where, 
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<== Cy = aa (D.8) 


Page letiing xX = r/a: 


Ww WwW 
O O 


2 oat oo pee lee eee = ’ = 2 ’ 
aoe = U + WoW = aa coum ee 2x) Cy = (2x 3X )C. 
ee Cee 2. 0 
u we 6 
~ _- = —_ iz J ete es ‘ 
2) a oa [ (1 x) Cy + (x-x )co] (D.10) 
and 
ee 3 7 
Sy 26, a a te - #x)C) + (3x - $x )Cy 
oS ee | 
+x (7 Ti* ea | (Died) 
Ue 
: 2 ae ay ’ ee ’ 
Eg + re. = aa [ (4 2x)c, + (2x +X )Co 
ree 2 Ces : 
mente < ) | a 


These expressions may now be substituted into the equations for 
moment and extensional force resultants, (2.14) through (2.17). 


Defining: 


Soe Gee, 
A. = 1 1 Sey 
5 O y Se 2 Wo 2,28 12 2,2 0 
ey * (2X ~ Bx) 5meD + GT Te)“ 
= ae a 30 2 
De a. il uf 7 
W W 2 W 
y= 1 _O ~ Sx2yor 28 4 X_ (28 _ 12,2)2_0 
E, = (4 2x)cy ot (2x Pe aa! > aT T1* ) °p 


oF 





we may write: 


« 
Meee A ei 2+1/n o, 1+1/n li j en 
ee Sk nfl (37) Isa IB, + 2A! 
_ n 1 2a lyn 1 _ 4] 7 1 
Tantiyaz!Bs + 25/5 * gat! Bs ~ 24s 
24+1 Pe 
+ Tosh aypaelee = 2A |. [D4 ( = x7] (D.13) 
Ss 
a 
Pe n 2+1/n o,l+l/n 1 bie eli/Ae 
eee 3K) nel” (37) A. (|B, ‘ 2A.| 
WwW W 
-|[Bo - 4al[et'/] - f(a - axes po + (2x - 3x?)e3 (9 
W 
2,28 _ 12.2)2 _o 
+ (TT ix) ™ J (Dai4) 
cmt lyn, Ovltl/n . 1 ae |) Ueda 
cee oe nen | (7) 2D. JE, + 3D.| 
7 n 1 27 it 4 it 7a 
(2n+1)D- JE. e EDA _ 2D. JE. De 
n es Ze ih r Zo a 1s 2 
2 (2n+1)D~ JE. oD |. TT Th see (D. 15) 
et A ene 1) en 2+1/n he 1l+l/n 1 1 1+1/n 
Uy ae SaTR amp (32) Do []E, + 2D] 
W W 
as 1+1/n fe ’ _O ice & 1 _O 
JE. 3D. | } [(1 x)C, 5 + (X-xX )cs - J 
(D.16) 


0 








The internal energy dissipation rate is: 


e _ 2 1 ® e e 
vVe= Zia f (M akg oo Ne, 


~ Ny&) XaxX 
giving with (D.13) through (D.16): 


W 
pe, a owl ry on 2+1/n 0, 1+1/n 
Vests” New aol Le (xz) : 


S 


fq é6w rdr 


(D.17) 


(D.18) 
where qT. 1s the integral of the bracketed geometric terms in 
the above expression. 


From the principle of virtual velocities 


(D.19) 
l/n ,l1 ,1+l1/n eae Tet Ln ee 
we (ee) 2ma* (35) h I. 6w = 
del rdr = anqée. 2e a2. (D+ 20) 
O a ia o 22 ; 
SO: 
q2t2/n a 
Yo 7 aif Cot ee Wie) 
3K S 
a. avetly me, ska, l/n it n 
5 S 
Se a,2nt+l 3 1 n 
Kaq _ am J 
| 5 S 


og 
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Illustration of Technique for Finding the Minimum Energy Integral 
As explained in Chapter II, the quantities Cy and C4 are 


determined by requiring that the internal energy dissipation 
rate (2.37) be minimized. Since the integration required to 
evaluate (2.37) must be done numerically, this minimization 


must be accomplished through a systematic search. The 
W 
Minimization must be done for each case (—, Mi OLetnterest- 


but, by considering a systematically varied series of cases, 
a set of curves such as those in Figures 2.2 and 2.4 can be 
obtained. 


The procedure used for minimization can be illustrated 
Ww 
by considering a typical case, n = 5.0 and = = 1.0 for 


clamped edges. The computer program listed in Appendix H, 
when provided an initial value for c) and a value for ci, 


will compute the integral, Tar for the given Cy and C5 values 


/ at regular intervals of 0.3. In the accompanying plots, the 


quantity ¢ rather than the integral I. is plotted so that $6 

may be transferred directly to Figure 2.2. As @ = £(=-), the 
E 

maximum value of ¢ corresponds to the minimum value of I: 


As seen in the plots, Cy is varied until a peak is found for 


a particular value of c Peaks are then found for other 


5° 
values of Cy until a maximum is found with respect to Ca. The 


resulting envelope of curves gives a maximum value of ¢ which 


gives a point on the plot of Figure 2.2. The procedure is 
W 


then repeated for the other values of = and n. 
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There is no mathematical assurance that the value of $6 
found as above represents a unique and absolute maximum value. 
The initial area of search was that suggested for the analytical 
solution of the n = 1.0 case; and a very broad parameter search 
with respect to Cy and C5 gave no evidence of any higher 
maxima. The agreement between the experimental results and 
the theoretical results suggests that the values of 4 obtained 


are, in fact, the prover ones. 


OZ 





APPENDIX F 


Tabulation of Data from Plate Experiments 


The following tables list the data taken during the plate 
experiments discussed in Chapter I.A. All tests were taken at 
room temperature which generally was 75°F +3°F with some points 
outside these limits. It is, of course, a well known fact that 
creep behavior in metals is affected by temperature. Within 
the relatively narrow temperature band here considered, however, 
such effects may be expected to be rather small and definitely 
secondary to the influence of pressure and geometry effects 
which were to be investigated in the plate tests. Examination 
of the data reveals no appreciable systematic influence of 
temperature ms these tests. All times are in hours and minutes, 
all deflections are in inches, all pressures in pounds per 


square inch, and all temperatures in °F. 
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feo 


8 LB. PLATE Sse 
Time Temp. 
Hrs-Mins r=0 r=+" r=1" r=14" r=2" r=2+" oF 
00 0052 IPD: 
one 53 .0049 .0033 .0021 .0005 .0001 
03 54 49 334 PA 05 013 
15 57 53 36 De, 06 014 
32 584 55 or 24 06¢ O14 qe 
50 60 57 384 24 07 02 
1-00 61 57 39 24 OF 02 72 
1-30 63 59 49 25 074 02 
2-00 64 61 414 26 09 02 72 
2-42 66 63 43 27 09 02 
3-00 67 64 43 274 09 02 72 
4-45 7. 67 61 29 10 02 
5-20 724 68 47 292 10 02+ 72 
6-00 73 684 48 30 10 024 
6-30 74 70 484 31 10 024 72 
25-00 92 88 62 39 14 07 74 
47-00 0110 0105 75 47 17 08 76 
52-00 14 08 78 48 18 084 76 
73-00 23 20 86 53 20 10 76 
100-30 30 264 90 554 21 10 76 
125-40 34 30 93 57 213 Tal 76 
TEST 2 
8 LB. PLATE 250 josul 
Time r=0 r=4" r=1]" r=]4" r=2" r=2¢" Temp. 
: - 
00 .0052 
01 53 .0045 .0036 .0020 .0007 .0002 Ta 
03 54 45 37 21 07 02 
15 56 51 384 23 0384 02 
32 58 53 40 234 09 02 76 
50 59 53 41 24 09 02 
1-00 594 534 41 24 09 02 76 
3-50 65 62 45 28 Ld 02 76 
4-30 67 64 46 28 1134 022 76 
5-00 68 642 47 29 12 03 76 
5-32 69 65 48 29 12 03 76 
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Time 


Hrs~-Mins 


Time 


8 LB. PLATE 


) Ibshe 


r=$" 


PLATE 


r=$" | 


TEST 


TEST 4 


r=1" 


SO Zon. 


3 


veel 


-O181 


4 


r=1s" 


0145 


4.1 psi 


r=2" 


-0060 
63 
65 
67 
69 
fe" 
73 
74 
75 
77 
80 
81 
83 
85 
86 
87 


4.1 psi 


r=2" 


~0055 
58 





Jlohey Woes: 


8 LB. PLATE 6.1 psi 
Time | | Temp. 
Hrs-Mins r=0 r=5" r=)" t= r=2" r=2-" ioe 
00 .0824 79 
01 .0869 .0805 .0614 .0376 .0150 .0055 
03 96 28 32 87 55 57 
06 .0916 47 46 95 59 59 
15 44 75 66 .0408 64 62 
25 61 Q4 81 5 68 64 
31 68 .0902 86 21 70 64 79 
45 87 15 98 7 74 66 
1-00 99 28 .0708 34 ay 66 80 
1-30 O20 48 23 44 80 68 
2-00 35 64 34 52 84 70 80 
3-00 59 89 53 64 90 Te 
3-36 72 1000 62 69 93 73 81 
4-45 93 20 76 79 98 75 
5-15 1107 2a 86 83 0200 82 81 
6-00 19 40 93 88 03 82 
TEST 6 
8 LB. PLATE 6.1 psi 
Time r=0 3 = r=}" r=)" r=]}" r=7) r=2-" Temp. 
oF 
00 .0840 
ont HOscwemmnoe42en.0631 9.0379  #.0156 .0073 76% 
03 .0908 65 49 90 61 75 
07 33 86 66 .0400 65 78 
15 55 .0910 83 11 70 80 
34 85 36 .0707 24 76 85 76+ 
48 98 49 Ly: 30 79 86 
1-00 . 1006 57 2D, 34 80 87 76+ 
1-25 21 70 35 4l 83 89 
1-50 33 81 44 47 87 92 qi) 
2-00 38 87 48 49 88 93 
3-15 61 .1009 68 59 92 95 78 
4-20 78 23 81 66 96 97 
5-50 98 43 98 76 PO Onn O00 78 
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Time 


Hrs~Mins 


Time 


8 LB. PLATE 


=0 = 3" 
1270 
Peis 25) 
39 76 
65 98 
Ger ncto2) 
1429 Sy 7! 
54 84 
75 .1403 
91 DI 
1524 49 
34 58 
52 76 
8 LB. PLATE 
r=0 r=+" 
1252 
1301 1247 
28 74 
43 88 
89 1333 
1413 54 
a 70 
45 83 
56 93 
65 1403 
89 23 
1504 40 
ly, a1 
24 57 
36 68 
50 33 
53 86 


ol 


r=]" 


O07 


TEST 7 


a oe mL 
r=l135 


-0613 


- 0586 


8.01 psi 


r=2" 


0244 


50 
55 
60 
70 
76 
81 
86 
93 
96 


20°30 


8.01 psi 


r=2" 


sz 


43 
47 
57 
62 
66 
69 
72 
i> 
80 
a5 
88 
89 
ee 
96 
2 


TEST 9 


6 LB. PLATE DASE) soul 
Time Temp. 
Hrs-Mins r=0 r=}" r=)" r=15" r=2" r=2¢" oF 
00 0099 74 
01 .0100 .0091 .0067 .0025 .0012 .0001 
09 LOS Seles 744 31 132 01% 
15 15 05 Val 32 14 012% 
21 18 074 78+ 33 142 02 
30 20 ta 80 343 15 02 
45 De 134 822 36 16 03 74 
1-00 252 Le 84 37% 16$ 03 
1-15 274 iby) 86 38 162 03 
1-30 294 192 874 382 1g 033 
2-00 322 232 90 402 18 05 74 
3-17 37 30 934 43 19 05 
4-00 382 312 95 44} 19 O54 74 
4~30 394 32 962 443 19 054 
4-45 392 a1, 964 443 19 05s 
5-00 41+ 33+ 97 454 19¢ 054 74 
5-15 41+ 334 97% 454 194 054 
5-30 414 332 9732 454 194 054 
6-15 412 362 99 47 20 054 74 
TEST 10 
6 LB. PLATE 230 Upsi 
Time = eet ei r=14" =o) ae a Temp. 
oF 
00 .0097 74 
01 -0100 .0089 .,0068 -0026 POCO. OOOO 
09 09 .0100 ae 31 12 00 74 
15 12 03 774 32 fe? 00z 
Di 14 05 79 33 13 Ol 
30 164 07 god 34 13 01 
46 194 11 got 36 132 01 74 
1-00 214 13 84 37 14 01 
1-15 23 134 85 38 147 01 
1-30 244 142 86 38 14% 01 
2-00 274 18 884 392 15 02 74 
215 33 23 924 42 16 o2t 
4-00 344 252 94 434 i 024 
4-30 35 Da 944 44 17 02% 
5-00 35% 28 95% 44% 17 022 
5-15 354 28 95% 44% le 022 
5-30 354 282 96 45 17 03 
6-15 354 29% 97 45% 17 03 





Time 
Hrs-Mins 


Time 
Hrs-Mins 


6 LBS Ee LATE 


6 LB. PLATE 


iow 


~0265 


TE See 2 


~0247 


UO 
86 
59 


Svour 
98 
01 
02 
05 
06 
10 
ys 
18 
19 
21 
23 
24 
25 
26 
264 


76 


76 


Time 


Hrs-Mins 


Time 


TEST 13 


6 LB. PLATE 6. | psi 
Temp 
va 8. r=4" cs) Oy r=14$" 2 r=2¢" or 
0889 76 
0927 0852 0650 0441 0175 0071 
51 qa 68 53 80 73 
68 90 80 62 84 75 
85 .0905 94 al 88 78 
94 15 .0701 76 90 79 
1017 38 20 90 96 81 76 
29 49 28 96 98 83 
33 54 33 99 0200 83 
45 64 42 0505 03 85 76 
57 78 53 ai 05 87 
66 84 59 15 06 89 77 
83 0100 71 24 AA: 91 
94 11 80 31 ae 92 77 
1104 22 89 37 15 93 774 
26 43 0807 50 Dil; 95 
33 51 14 54 23 96 77s 
45 60 23 60 26 0100 
TEST 14 
6 LB. PLATE 6.1 psi 
r=0 Tse r=]" r=1;" r=2" Sa aey Temp 
oF 
0878 
0936 0865 0643 .0438 .0152 .0078 80 
62 88 65 51 57 82 
87 0915 85 65 62 86 
1016 42 0706 79 68 97 80 
Be 56 18 88 yal 98 
44 68 28 94 74 .0100 80 
60 83 40 .0503 78 03 
66 88 45 06 79 04 
79 98 55 1s 81 05 80 
Ooms 66 Da) 84 07 
97 a7 70 24 84 07 
1106 26 i 30 86 08 80 
18 37 87 36 89 09 
30 48 98 44 92 11 
36 56 0802 48 94 12 80 
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6 LB. PLATE 8.0 psi 
Time 
Hrs-Mins r=0 r=4" ssh r=14$" r=2" 
00 -1308 
Ol ~1354 1265 1009 0679 20255 
03 78 88 28 93 61 
06 98 1308 44 0705 66 
1} 1427 36 a2 Zo 74 
24 57 54 84 36 80 
30 63 61 92 42 83 
45 We 78 1108 54 88 
1-00 87 91 19 63 92 
1-15 97 1403 29 70 96 
1-36 fe Bea 19, a4 40 79 0300 
2-00 24 Pye Sa 88 05 
2745 44 Se 69 0803 Jak 
3-00 50 59 74 07 12 
4-20 74 81 97 Zo 20 
5-04 87 95 1206 33 24 
5-30 93 .1500 12 Bis 26 
6-00 1600 06 18 43 28 
TEST 16 
6 LB. PLATE 8.0 psi 
Time r=0 r=}" r=1" r=]4" ioe 
00 1295 
Ol ~Il342. 71246 1009 -0636 ~0265 
03 68 78 30 53 72 
06 88 97 47 65 vay! 
15 lea sey si 322 73 83 86 
24 34 44 88 oF 91 
45 65 68 AS 0716 0300 
L700 78 85 25 27 05 
1-16 90 98 37 35 09 
1-30 99 - 1406 44 47 dN 
2-00 1517 19 60 53 18 
2-30 ai 36 71 63 23 
3-04 44 47 83 Ve 27 
3-30 53 59 92 80 30 


Meo 15 


iL 


76 





Time 


Hrs-Mins 


Time 


4 LB. PLATE 
r=0 3 _ r=" 
0110 
0112 0106 
Ne) 08 
18 13 
AAS) iG 
27 21 
33 26 
39 a2 
45 38 
sys) 46 
56 48 
57 50 
4 LB. PLATE 
r=0 r=$" 
0127 
0131 OF22 
S15) 25 
35 2 
59 eal 
42 Sy) 
45 oy 
47 oo 
49 42 
AG 42 


ES 17 


eZ 


jos | r=ls 
0080 .0050 
82 50 
84 Sy 
88 53 
91 56 
94 58 
98 59 
0103 63 
08 66 
10 67 
ae 68 
IMR SIE ALS: 
je= lL fess lle A 
-0080 .0055 
81 56 
83 57 
86 59 
87 60 
89 61 
89 62 
91 63 
91 63 


-0018 
18 
20 
20 
21 
22 
23 
25 
26 
27 
OR 


-0020 
Psd 
23 
22 
23 
24 
24 
25 
25 


Temp 
ain (Op 
ge 
03 
04 
04 74 
04 
04 LS 
05 
05 74 
06 
06 74 
06 
Temp 
a 
2S 
19 
20 
20 v3 
21: 
22 Vee: 
22 
22 
a 73 





WES l oo 


4 LB. PLATE 1.840 psi 
Time Temp 
Hrs-Mins xr=0 r=}" r=" 2 ee oe oP 
00 .0845 76% 
O1 .0871 .0820 .0653 .0452 nou Oe 0045 
03 83 32 64 55 73 46 
06 92 Al 72 61 75 47 
15 0906 55 83 68 79 48 
21 11 60 88 71 81 48 
31 18 68 93 76 83 49 764 
45 24 74 98 80 84 49 
1-00 29 80 .0702 83 86 50 764 
1-31 34 86 08 88 87 50 
2-00 42 92 13 92 89 51 764 
2-30 47 98 18 95 91 51 
3-00 51 .0901 22 98 92 52 a 
TEST 20 
4 LB. PLATE 1.840 psi 
Time Te) ae fees pedigy 2) ce weer=2.° > Temp 
oF 
00 .0820 79 
01 .0848 .0783 .0638 .0438 .0162 .0053 
03 61 98 49 46 65 55 
06 70 .0804 56 50 68 56 
12 79 14 63 56 70 57 79 
18 85 24 68 61 VY 58 
26 91 31 qo 65 74 58 794 
30 94 33 75 67 75 58 
1-00 .0907 46 85 74 79 61 794 
1-15 1: 49 89 76 80 63 
1-30 16 56 92 79 81 63 79 
2-00 20 58 98 84 83 64 
2-30 28 67 .0702 86 85 65 
3-30 35 76 09 92 88 66 79 
4-19 44 84 14 96 91 68 
4-30 44 86 MS 97 92 68 79 
5-30 50 92 19 .0500 94 68 
6-00 52 94 Ball 02 95 68 79 
6-33 55 96 22 03 95 68 


13 





Tame 


Hrs-Mins 


Time 


4 LB. 


PLATE 


TEST 21 


4.1 psi 
Temp 
sa a eae Tee oe sa 
794 
1333 0998 0377 0159 
1G ealioale 83 61 
57 19 88 64 
72 31 95 67 
77 34 97 69 794 
87 44 .0402 70 
90 47 04 gpl 80 
1403 57 10 74 
06 61 ll 75 80 
16 68 16 79 
23 75 19 80 80 
30 81 24 82 
32 83 25 83 80 
a7 88 28 87 794 
41 91 Bul 89 
45 94 32 90 79- 
TEST 22 
4.1 psi 
me Eels Ge DEON po Temp 
oF 
81 
1313 .0973 .0435 0136 
29 85 42 40 
40 95 48 42 
56 .1011 56 46 
65 19 63 48 
69 23 65 49 
79 32 eal 51 
83 35 73 52 81 
92 43 78 54 
99 50 82 55 81 
mal: 59 86 59 
17 66 90 60 81 
20 69 95 62 
23 gal 96 62 804 
28 75 98 63 80 
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Lol 23 


4 LB. PLATE 6.1 psi 
TIME Temp 
Hrs-Mins r=0 r=i" r=l" r=1)" r=2" r=2-" 219 
00 CIES: 76 
Ol Poo. toz0me. 1554 .1138 -0488 .0243 
03 39 35 70 53 96 46 
06 51 47 83 63 0503 Sal 
15 66 66 95 oh 12 54 
209 WES) 73° .1601 85 18 55 76 
3] BZ 85 07 onl DA, 57 
45 85 89 09 Os 1, § 58 
= 0 0 89 94 13 94 . 26 60 76 
15 wes! 94 ily, 2) 27 61 
t= 33 95 PO L Sets ZO? Zo 61 
200 oo ON 23 05 30 62 76 
3-00 - 2009 07 20, ON, 34 63 76% 
4-00 14 08 30 12 bbe, 64 Ti 
5-00 ig 18 39 17 38 fall 
a3 20 19 a9 18 a9 71 774+ 
TEST 24 
4 LB. PLATE 6.1 psi 
Temp 
Time r=0 r=" r=)" r=)+" r=2" r=2¢" Oe 
00 Belles iss 84 
Ol PEO eeeo7 Seb 5S3 .1187 pU45L 20212 
03 - 2000 Wie lOOs «1203 62 18 
06 16 .1908 20 14 Al 22 
08 24 14 a7 20 76 24 
15 45 Si) 43 33 84 28 8 4 
26 64 60 58 47 o5 a2 
30 68 62 62 49 96 Se) 84 
45 81 TZ a3 58 -0502 36 
= 0'0 90 81 Bel: 13 07 38 84 
Nes S10) 2103 93 28) 73 LAS: Al 
2-00 14 .2004 1702 80 20 44 84 
2=54 24 10 10 85 IES) 48 
3-00 29 19 als 89 27 49 84 
2-350 35 20 20 94 30 50 
4-02 40 29 PRS AE SHONO. Si 51 84 
4-30 45 33 Zo 03 34 52 
Dali0 49 39 34 06 36 54 
byob 18) 54 44 By 09 38 55 84 
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APPENDIX G 


epect rochemical Analysis Report 


The reported analySis was made by Mr. Walter W. Correia 
of the STRNAD Spectrographic Laboratory at M.I.T. The sample 
labels, "4", "6", and "8" are for the 4 lb. plate, 6 lb. plate, 
and 8 lb. plate respectively. The analysis samples were 
taken from tensile test specimens randomly selected from 
those actually used in the tensile experiments. From this 


data one may infer: 


4 lb. plate Pine) 99.67% 
& tboea plate Ph > 99.96% 
8 lb. plate Poh > 99.293 
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APPENDIX H 


Computer Programs 


I. Program for Clamped Edge Case 

The following listing is for the program used in the 
numerical evaluation of the energy dissipation rate integral 
for the clamped edge circular plate. The program consists of 
a main part and two subroutines. Subroutine SAMS evaluates 
the moment terms, Me. + Mak gs at a radial position x, where 
x = rfa. Simlarly Subroutine JOES evaluates the extensional 
force terms, Né, + Npegs at a position x. The main program 
uses five point Gauss quadrature to evaluate the total 
integral over x. The program is written in FORTRAN IV. 

The inputs to the program aren, Cy, cy,wi/h, and Qs. The 
input format is F10.5 which means each input quantity is 
alloted 10 spaces and if no decimal is inserted by the 
programmer the decimal is assumed to be to the left of the 


last five spaces in the field. The format for the invut card 


is as follows: 


Freld Quantity 
0 n 

i= 2.0 Cy 
21-30 C5 
31-40 wo/h 

ABS 3 0 chy 


The last data card must have a negative number in the field 


1-10. This causes the program to exit. For each data card 


Ps 





Enews program Comoutes and prints seven sets of output data 
where: 

' — 
Saad 


starting with the input value of Cy olus .1l1. All other 
parameters remain constant through all seven iterations. 


The output is in ten columns as follows: 





Column Quantity 
1 n 
i] 
2 Ch 
? 
3 Co 
4 wi/h 
5 ay 
6 | AL! the integral of 
Ne as ; 
Se ae + Mako 
7 An? the integral of 
NECe + Noe 5 
8 Rah ea | 
n C 
9 ® 
10 go” 
3-04 
where 9 = 1D rT 
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NYOLse 

X*¥CC(HALISEVIF0 (HWS )SeV)) = XWA 
C2VRCESH( (OP CHHXRNO THT V HCP 7) Dewl = Bal 
C(eVs5—-2V Ed) x83 KC (Tvsa3t vas) (Cv a808Z)/0°T)) = wl 
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CINS 


NaoOLsaS 

X*¥(C(CSNLISHV)+((SNH)S9V)) = XNA 
SDHe#(NVSDRENVAOSD)#(VD/S0O°T) = SN1 

(HMe Ze ( (O° SeHX) @X) D+ (HMe TDe(X—CPT)) = JD 
(NS/O°T+O0°T) ##V5D = NVSD 

(39)S98v = vw3Dd 

DI=eKD = 3D 
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(CVHOTE VHC (CCE CHHX) HOF GRMOFE)HRTIVHO°Z) = 
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II. Program for Simply Supported Case 


The following program for the simply supported case is 
exactly like that for the clamped case except that a different 
assumed velocity profile is used. 


The input data, in format F10.5, is as follows: 


Field Quantity 
1-10 n 
11-20 Cy 
21-30 Co 
31-40 Wwofh 


Here again the last data card must have a negative number in 
Field 1-10 to cause the program to exit. The output data is 
exactly as in the clamped edge program discussed in Part A 


except that the column for Oly is omitted and here 
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QUADRATURE 
See 2369 ZE6E9S1#DI14+2.23692653514*D012 
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APPENDIX I 


Computation of Total Deflection 


The discussion in Chapter II was centered about determining 
the deformation rate, Woo In practical problems it is 


frequently the total deflection, Wor at some time, t*, which is 


of interest. It has been shown that ce is a function of =: 


and also that the material creep constants, K and n, may 
change with time. To compute Wer therefore, a piecewise 
linearization technique can be employed in which the deformation 
process is broken into a sufficient number of linearized time 
segments to provide the accuracy desired. One such technique 
is discussed in this Appendix. 

In dividing the:deformation process into time segments 
some guidance can be obtained from the Shape of the tensile 


: W 


test curves and from the shape of the 6% versus —— curves, 


Figures 2.2 and 2.4. The time segments should be relatively 
short during primary creep and may be longer when deformation 


is proceeding more slowly. As the deformation process 
W W 

One: O 
h increases; and when a 


a significant change in ¢? a new linear segment should be 


proceeds, has changed enough to make 
introduced. 

Having divided the total time of interest, from t=0 to 
t=t*, into n segments one can proceed as follows: 

Let t=t. be the time at the beginning of the segment 


under consideration and t=t. the time at the end. From 


£ 


tensile test data determine the appropriate values for 


ia 


Bieweueconecoenscants, K and n, at the time tae in the 


center of the given time segment, where: 


{ee (cima) 





Using these values of K and n; and using the value of 
Ww 


= at the beginning of the time segment, compute we 

with the appropriate equation from Chapter II, (2.53) for 
clamped edge or membrane or (2.57) for simple support. 
Multiply this value of Wa by the total time length of 


the segment, t and add the result to the value of 


fg! 

wy at ta. This gives Wo at ty. The procedure is then 
repeated for the next segment, and the next, and so forth 
until t* is reached. 

In illustration, consider the case of Test 13, a 6 lb. 
plate under 6.1 osi pressure. Assume that it is desired to 
compute We at t=8 hours. The total time is broken into two 
segments, a relatively short initial segment from t=0 to t=2 
hours, and a longer segment, where the deformation is more 
gradual, from t=2 hours to t=t*=8 hours. Values of K and n 
are determined at the median times of the segments, which in 


Ese case are t= hours and t=5 hours. From the tensile test 


results in Chapter II the appropriate values are seen to be: 


t=] hour n=5.0 meant ce eas (222) 

c=5 hours n=5.0 K=1.7 x Nowecs (i oy 
In test 13 the initial deflection was .0880" giving: 

x) = .840. (I. 4) 

h 


WS 





If test data is not available the initial deformation can be 
computed using elastic or plastic analysis. Using (I.2) and 
(I.4) and the appropriate value of © from Figure 2.2 in 
equation (2.48) the value of Wi for the first segment is found 


to be .0090"/hr. Deflection at t=2 hours may now be determined: 


(Wo ee =se0ee0 + 2(.0090) = .1060 5) 
which gives: 

Wo 

Srna nee. = 1.01 (I.6) 


Using (1.3) and (1.6) with Figure 2.2 and equation (2.48) we 
find for the second segment that ie = .0020"/hr. Deflection 
at t=8 hours is: 

Wo) bag 
A plot of this two step linearization and of the measured 


= .1060 + 6(.0020) = .1180" (ie) 


experimental curve is shown in Figure I.1. Also shown is a 
three step linear approximation made with the technique 

eee. tbed above with periods 0 < € < 1 hour, I < t < 3 hours, 
and 3 < t < 8 hours. Still closer approximations to the 
actual curve could be made by dividing the time period into 
more segments. 

In the practical application of these results it is 
important that initial deflection of the plate be known. 
Computing this initial deflection, however, may be no simple 
matter. Finite deflection analysis techniques for plastic and 
elasto-plastic deformation of plates are usually quite complex. 
For materials in which creep is important there is often a 


high degree of sensitivity to the speed of application of load 


1eS:3 


[28] which makes even the selection of appropriate elastic 

and plastic material constants a rather difficult task. In 
fact, such elasticity concepts as yield stress and Young's 
modulus may become rather meaningless for some materials under 
conditions where creep is important. 

In cases where the initial deflection is expected to be 
quite small and elastic, one could calculate initial deflection 
using the bending solution of Timoshenko and Woinowsky- 
Krieger [29]. For larger elastic deformation the large 
deflection solution of Timoshenko and Woinowsky-Krieger [30] 
is available. The large deflection elasto-plastic deformation 
of circular plates is treated by Ohashi and Murakami [31] [32], 
and the large deflection of elastic-perfectly plastic circular 


plates has been analyzed by Crose and Ang [33]. 
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